LATTICE STATISTICS ON PLANE, CYLINDER AND TORUS 

ZUHE ZHANG 



Abstract. In this paper the relations among entropy constants of lattice statistics on 
plane, cylinder and torus are discussed. We show that for many statistical models, 
the entropy constant of a plane lattice is the same as the entropy constants of the 
corresponding cylindrical and toroidal lattices. As a result, we get the entropy constants 
of quadratic lattices on cylinder and torus. Furthermore, we consider three more complex 
lattices which can not be handled by a single transfer matrix. Using the concept of 
transfer multiplicity we extend results of quadratic lattice to these three lattices and 
obtain lower and upper bounds of the entropy constants. 



1. Introduction 

The lattice statistics in statistical physics has a long history. A typical problem is to 
count the ways of putting particles in a plane lattice sites so that no two share the same 
site or are in adjacent sites. Such problems are called the planar lattice gases models 
[H El El m [HI [I3]. Mathematicians formulated them by the enumeration of the (0, 1) 
matrices which describe the independent sets in a plane quadratic lattice graph (also 
called a planar grid graph) [H [HI [lEl [17] . In this paper, we use Gm,n (where m, n always 
denote positive integers) to denote a plane lattice graph, that is a finite part of a plane 
lattice, whose vertices are arranged in (m + 1) rows and {n + 1) columns. Given an 
independent set S of graph Gm.n, a portion of S that lies in a fixed column of Gm,n can 
be represented by an (m + l)-vector of O's and I's, where a 1 indicates the vertex is in S 
and a indicates the vertex is not in S. 

If Gm,n is a plane quadratic lattice graph, then any (m + l)-vector arising this way has 
no two consecutive I's. Clearly, a vertex subset S of Gm,n is an independent set only if 
all its n + 1 corresponding (m + l)-vectors represents independent sets. Let Pm denote 
the set of all (m + l)-vectors of O's and I's with no two consecutive I's. The cardinality 
of Pm is well known (and can be easily seen) to be F^+s, the Fibonacci number (starting 
from Fq = 0, Fi = 1.) We can construct an (m + 1) by (n + 1) matrix M to represent an 
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independent set S of Gm,n by the gluing procedure described below. First, take a vector 
from Pm such that it corresponds to the first column of S, and denote it as Vi. Then, 
to the right of vi we glue a vector V2 selected from Pm that corresponds to the second 
column of S. Then we glue a vector ^3 from Pm to the right of V2, and so on so forth. 
Continue this way until the {n + l)-th column is glued and then we obtain the (m + 1) 
by (n + 1) matrix M representing the independent set S. We can get all the possible 
matrices representing independent sets of Cm,™ by this procedure of gluing vectors from 
Pm- Note that two vectors of Pm can be glued together if and only if they have no I's in 
common position, i.e., they are orthogonal (their dot product equals zero.) In the above 
procedure, we glue columns from left to right. Similarly, we can have another procedure 
that glues rows from top to bottom. Note that for non-grid lattice graphs Gm,ni these 
two procedures may lead to different transfer patterns. But for all the lattices considered 
in this paper, the transfer patterns in the two procedures are the same. 

Fig{T] shows an independent set S in the plane quadratic lattice graph G^^^. The por- 
tions of S that lie in each of the columns can be represented by the respective 5-vectors 
(0,0, 1,0, 1),(0,1,0,0,0),(0,0, 1,0,0), (1,0,0, 1,0), (0,0,0, 0,1),(0, 1,0, 1,0). 



Figure 1. 

From a plane quadratic lattice graph Gm,n {n > 1), identifying its m edges on the 
left with its m edges on the right, correspondingly, we get a cylindrical quadratic lattice 
graph Hn,m (Note that here the graph can be seen as drown on a vertical cylinder.) Let 
Gn denote the set of n-vectors of O's and I's with the property that no two consecutive 
I's occur in cyclic order. Then, similar to the discussion on Gm,n, we can see that any 
independent set S of Hn,m can be represented by an (m + 1) by n matrix M whose rows 
are from the set C„, and that all the representing matrices M can be obtained by the row 
gluing procedure (from top to bottom) similar to the column gluing procedure (from left 
to right) for Gm,n- It is not difficult to see that the cardinality of the set Gn is F„_i + 
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Similarly, by identifying the top and bottom boundary cycles of a cylindrical quadratic 
lattice graph we get a toroidal quadratic lattice graph Sn,m that can be seen as 

drawn on a torus. 

In general, for any plane lattice graph Gm^n, the cylindrical lattice graph and the toridal 
lattice graph obtained by the identifications as above will also be denoted Hn,m and Sn,m- 

For a given plane quadratic lattice graph Gm,n, the transfer matrix is an x Fm+3 
matrix of O's and I's, defined as follows. The rows and columns of are indexed by 
vectors of P^, and the entry of in position (a,/?) is 1 if the vectors a,f3 are orthogonal, 
and is otherwise. Note that the matrix depends only on m. 

For example when m — 3, the possible column vectors of are (0, 0, 0, 0), (1, 0, 0, 0), 
(0,1,0,0), (0,0,1,0), (0,0,0,1), (1,0,1,0), (1,0,0,1), (0,1,0,1). If we index the rows and 
columns in this order, then the transfer matrix of Gm,n is 
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Similarly, for cylindrical quadratic lattice graph H^^jn, its transfer matrix T„ is an 
{Fn-i + Fn+i) X {Fn-i + -^n+i) Symmetric matrix of O's and I's, defined as follows. The 
rows and columns of are indexed by vectors of C„, and the entry of T„ in position 
{a, (3) is 1 if the vectors a, (3 are orthogonal, and is otherwise. For n — A the possible 
row vectors of are (0,0,0,0), (1,0,0,0), (0,1,0,0), (0,0,1,0), (0,0,0,1), (1,0,1,0), 
(0, 1, 0, 1). If we index the row and column in this order, then the transfer matrix of i^n.m 
is 
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For a plane quadratic lattice graph Gm,n, by adding two crossed diagonals in each square 
inner face we get a crossed quadratic lattice graph, which will be denoted as Gm,n- The 
transfer matrix of Grn,n can be defined in a similar way. We say that Gm,n and Gm,n are 
lattice graphs with transfer multiplicity one, since computing the number of independent 
sets for each of them only needs to employ one transfer matrix. This may not hold 
for other lattices. For example, in [19], for an aztec diamond we need to introduce two 
transfer matrices, each of which is the transpose of the other. So the transfer multiplicity 
of an aztec diamond is two. There are also lattices with transfer multiplicity three; one 
example is the 8.8.4 lattice that will be discussed later in section 4. 

The entropy constant of a plane lattice is defined hy rj = lim /(m, n)^/*^*-™'"-' where 
/(m, n) denotes the number of independent sets of Gm,n and fc(m, n) denotes the number 
of vertices of Gm,n- The entropy constants of cylindrical and toridal lattices can be defined 
similarly. As in [1], we have /(m,n) = ^ T"^ = 1 ■ T"l for G^.n and i?n,m- Clearly, 

Trace(T")= ^ T:^u- 

In |1] , Calkin and Wilf proved the existence of the entropy constant of plane quadratic 
lattice and established its upper and lower bounds. Two natural problems are to consider 
the entropy constants for lattices on cylinder or torus. Note that the method of Calkin and 
Wilf 's is valid for the lattices with the same symmetric transfer matrices in both horizontal 
and vertical directions. In this paper, we will investigate the relationship among the three 
cases (plane, cylinder and torus) for several lattices. We will show that for each type of 
four lattices, the entropy constant is the same no matter the lattice is on plane, cylinder 
or torus. Upper and lower bounds of the entropy constant will be established for the 
lattices discussed in this paper. 



lattice statistics on plane, cylinder and torus 
2. Lattices with transfer multiplicity one 
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The entropy constant of the plane quadratic lattice was already discussed in [4]. Now 
we consider the crossed quadratic lattice, which is obtained from the plane quadratic lat- 
tice by adding two crossed diagonals to each square inner face. Figl2] shows the crossed 
quadratic lattice graph 64^5 that is a part of the crossed quadratic lattice, where an 
independent set S is indicated by small circles. The portions of S that lie in each 
of the columns can be represented by 5-vectors (0,1,0,0,0), (0,0,0,0,1), (1,0,1,0,0), 
(0,0,0,0,1), (1,0,1,0,0), (0,0,0,0,1). 




Figure 2. 

Note that the plane crossed lattice graphs are no longer planar graphs. However, it 
is easy to see that the transfer matrix of the plane crossed lattice graph Grn,n is an 
Fm+3 X -Fm+3 matrix of O's and I's. When m = 3, the possible column vectors of Pm are 
(0,0,0,0), (1,0,0,0), (0,1,0,0), (0,0,1,0), (0,0,0,1), (1,0,1,0), (1,0,0,1), (0,1,0,1). If 
we index the rows and columns in this order, then the transfer matrix of Gm,n is 

l\ 







Clearly, the plane crossed lattice is a non-planar lattice with transfer multiplicity one. 

For a given Gm,n, by identifying its left edges with the right edges correspondingly, 
we get the cylindrical crossed lattice graph Hn,m- The transfer matrix of Hn,m is an 
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{Fn-i + Fn+i) X {Fn-i + i^n+i) Symmetric matrix of O's and I's. When n — A the possible 
row vectors of its C„ are (0, 0, 0, 0), (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1), (1, 0, 1, 0), 
(0, 1, 0, 1). If we index the rows and columns in this order, then the transfer matrix of 

Hn,m is 
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In the following lemma, we will establish a relation between the transfer matrix of 
Gm,n (cylindrical lattice Hn,m) and the number of independent sets of the corresponding 
cylinder lattice graph Hn,m (toroidal lattice graph Sn,m)- Then we will prove our main 
theorem and provide numerical upper and lower bounds of the entropy constants for the 
quadratic and crossed lattices on cylinder and torus. 

Lemma 2.1. For each plane lattice graph Gm,n with transfer matrix T^, transfer multi- 
plicity one and a positive integer n, the trace of is equal to the number of independent 
sets of the corresponding cylinder lattice graph Hn^rn- For each cylinder lattice graph Hn^ra 
with transfer matrix T^, transfer multiplicity one and a positive integer m, the trace ofT^ 
is equal to the number of independent sets of the corresponding torus lattice graph Sn^^. 

Proof. Recall that Hn^m can be obtained by identifying the left column and the right 
column of Gm,n- Thus there is a bijection between the independent sets of and the 
independent sets of Gm,n whose left and right column vector are the same. And the latter 
is the trace of T^. 

Similarly, Sn,m can be obtained by identifying the top row and the bottom row of i?n,m- 
Thus there is a bijection between the independent sets of Sn,m and the independent sets 
of i?n,m whose Corresponding top and bottom row vectors are the same. And the latter 
is the trace of T^. □ 

Theorem 2.2. For a lattice with transfer multiplicity one, if in both directions the transfer 
matrices are the same real symmetric matrix, then its entropy constants on plane, cylinder 
and torus are the same. 
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Proof. Let T be the transfer matrix of columns of Gm,n with the characteristic polynomial 
/(x), and let 

where g{m) equals the number of independent sets on the first column of Gm,n- Particu- 
larly, for the quadratic lattice and the crossed lattice, g{m) = Fm+s- 

By Hamilton- Cayley Theorem we have 

/(T) = aol + aiT + aaT^ + ■ ■ ■ + a,(„)r3(-) = 0. 
Put /o = 1 and /„+i = T/„, then 

fn = Tn. 

Let 

'^m,n 1 ■ fn 1 ' X' 1, 

where 'Um.n is the sum of all entries of the matrix T", namely the number of independent 

sets of Gm,n- 

Thus 

1 ■ f{T)l = CloUmfl + aiUm,l + a2'Wm,2 + " " " + ag(m)Um,g(m) = 0. 

In general Um,n satisfies the following recurrence relation: 

-|- ■ ■ ■ -|- CLg(^m)'^m,n+g{m) 0- 

By a well known theorem on difference equations (see, for example, [H]), the charac- 
teristic polynomial of this linear recurrence relation can be written as 



f(x) = {x- XiY^{x - X2Y^ ■ ■ ■ (s - A,)^- where Ai < A2 < • • ■ < A, 
Note that A^ depends on m. Then 
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Um,n = Y.Pi{n)Xl 
1=1 

for all n where Pi{n) is a polynomial with degree at most — 1 in n. The coefficients 
Pi{n) of the polynomial Um,n are determined by the initial values 

^m,0^ ^m,l ) ' ' ' ) ^m,g{m) — l ■ 

Note that the first row and the first column of T are both vectors of all I's. Thus the 
matrix T is non-negative, irreducible (its corresponding digraph is strongly connected) 
and prime. Hence the spectral radius of T is a simple positive real eigenvalue with 
magnitude greater than any other eigenvalues' (see [15]) and Ps(n) is a positive constant. 
It is not difficult to see that 



lim («„,„) = A^. 

n— ►oo 

Since in both directions the real and symmetric transfer matrices are the same, by the 
same method in [1], one can show that the following double limit exists: 

Tji = lim (um,n)^^™ = lim A^. 

m,,n— ►oo ' m— soo 

Now by Lemma 2.1 we can see that the trace of the n-th power of the transfer matrix 
of Gm,n is equal to Vm,n^ the number of independent sets of -ffri,m- 

Since the trace of is the sum of its eigenvalue of T" and each one is an n-th power 
of an eigenvalue of T, we see that 

if 

s 

Vm,i = TraceiT) = X^CiAi, 

i=l 

then 

Vm,n = TraceiT^) = Ee^r, 

i=l 

where Cj is the multiplicity of eigenvalue A,. 
Since A^ is the simple largest eigenvalue, then 

lim (wm.n)^/" = K- 
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So, the entropy constant of Hn,m 

T]2 = lim = lim xV"" = r],, 

m,n— »oo m— »oo 

That is, the entropy constant of Gm,n and Hn,m are the same. Similarly, using the 
second conclusion of Lemma 2.1 we can prove that entropy constants of Hn,m and Sn,m 
are the same. This completes the proof. □ 

In |1] Calkin and Wilf already obtained a good estimate for the entropy constant of the 
plane quadratic lattice. Now, by Theorem 2.2 we immediately see that the entropy con- 
stants of the quadratic lattice on plane, cylinder and torus are all between 1.503047782... 
and 1.5035148.... 

Since the transfer matrix of crossed quadratic lattice is symmetric, the same approach in 
[4| and the proof of theorem 2.1 can be processed here. Thus we can use the same method 
in [Ij to get the upper and lower bound of the entropy constant of crossed quadratic lattice. 
The lower bound of lim f{m, where A's are the largest eigenvalues 

m,n—>oo 2g 

of corresponding T's and the upper bound of lim f{'m,nY^"^"' is (i^2fc)^^^'^ where ^'s are 

r7i,n— >oo 

the largest eigenvalues of corresponding B's. Let p = 4, q = 4 and k = 6, we have 
1.342542258... < t] < 1.342652572.... 

3. Lattices with transfer multiplicity two 

In this section we consider the lattices with transfer multiplicity two, i.e., the lattices 
for which we need two transfer matrices to compute the number of independent sets of 
the lattice graphs. One of such lattices is inspired by the famous aztec diamonds. The 
study of enumeration of perfect matchings, spanning trees and independent sets of an 
aztec diamond can be found in [5l [6l [12] and the references cited therein. 

Let Li be the path with i vertices 1, 2,- ■ ■ ,i. The tensor product of two paths ® is 
the graph on m X 77, vertices {{x,y) : 1 < x < m,! < y < n}, with {x,y) adjacent to {x\y) 
if and only if |a: — = \y — y'\ = 1. This graph consists of two connected components, 
the one with the vertices {(x, y)\x + y is odd}, denoted 0{L2m+i ® -^^2m+i), is called the 
aztec diamond of order m. More generally 0(L2m+i ® -^2n+i) is called the generalized 
aztec diamond of order m x n introduced by the author in [19]. An independent set of 
a generalized aztec diamond can be represented by an ordered list of column vectors. 
FiglS] shows an example. For 0{Lq (g) Lg), the given S (the big dots) can be represented 
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by the ordered list of 9 column vectors: (0, 1, 0, 0), (1, 0, 0, 1, 1), (0, 1, 0, 0), (0, 0, 0, 1, 0), 
(1,1,0,0), (0,0,0,0,0), (0,0,0,1), (0,1,1,0,0), (0,0,0,1), (1,1,1,0,0). 




Figure 3. 

By identifying the top row and the bottom row of a generalized aztec diamond, we get a 
cylindrical generalized aztec diamond. And, the toroidal generalized aztec diamonds can 
be obtained by identifying the left cycle and the right cycle of a cylindrical generalized 
aztec diamond. 

Now we consider the transfer matrix of 0(L2m+i®-^2n+i)- It is clear that the generation 
of each independent set of 0{L2m+i <S> L2n+i) involves 2n assembling steps after the first 
column Vi is established (For simplicity, the assembling of the {i + l)-th column f j+i to 
the ith column Vi is called step i.) Step one is to assemble V2 to the right side of vi. 
The transfer matrix representing step one, denoted T^j, can be constructed as follows. 
Let Rm be the set of all possible vectors Vi. Obviously, Rm consists of all m-dimensional 
vectors of O's and I's, so it has 2™ vectors. Similarly, the set of all possible V2 is the 
set Rm+i of all (m + l)-dimensional vectors of O's and I's, and Rm+i has 2™+^ vectors. 
Then the transfer matrix = {Tvi,v2) is a 2*" x 2™+^ matrix whose rows are indexed by 
vectors of Rm and columns are indexed by vectors of Rm+i, where T^^^^j = 1 if f i and V2 
represent possible consecutive pair of columns in an independent set of 0{L2m+i® L2n+i), 
and T^,^^t,2 = otherwise. Similarly, the transfer matrix for step two is a 2"^'^^ x 2™ matrix 
which is the transpose of Tm^- It is easily seen that Tmi is the transfer matrix for 
every step i where i is odd and is the transfer matrix for every step i where i is even. 
Thus, if we take the transfer matrix of generalized aztec diamond to be Tm = TmiTm2 
then it is a 2*" x 2™ symmetric matrix. Furthermore, the transfer multiplicity of can 
be considered as one and the results of Theorem 2.2 hold. 

When m = 3, if we index the rows and columns in increasing order in binary numbers, 
then the transfer matrix of Gm,n is an 8 x 8 matrix given by the product of T^^ and T^^ . 
Here 
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Let Rn denote all n-dimensional vectors of O's and I's. The cylindrical generalized 
aztec diamond, obtained by identifying the left and right columns of the generalized aztec 
diamond 0{L2m+i <8) L2n+i), can also be seen as obtained by beginning with some vector 
in Rn, gluing to the top a new one in Rn such that the vertices represented by the I's 
in these two vectors are not adjacent until 2m + 1 vectors are glued. It is clear that for 
each i, step i and step i + 2 can be represented by the same transfer matrix. Let 
{Bn2, resp) denote the transfer matrix of the cylindrical generalized aztec diamond, which 
represents every step i for i odd (even, resp). Clearly, is a 2" x 2" matrix of O's and 
I's whose rows and columns are indexed by vectors of Rn, the entry of in position 
(a,/?) is 1 if the vectors represent possible consecutive pair of rows in an independent set 
of 0{L2m+i ® L2n+i) ou Cylinder, and is otherwise. It is no difficult to see that Bn2 is 
the transpose of B^- Thus we get the transfer matrix of cylindrical aztec diamond as 
Bn — Bn^Bn^, which is a 2^* X 2" symmetric matrix. Furthermore, the transfer multiplicity 
of Bn is one. 

Take n = 3 as an example. If we index the rows and columns in increasing order in 
binary numbers, then the transfer matrix of Hn,m is an 8 x 8 matrix given by the product 
of B^^ and B^^ , where 
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Since the generalized aztec diamond lattice has the same symmetric transfer matrices 
in both horizontal and vertical directions, the same approach in |1] can be taken here. So, 
the lower bound of lim /(m, tt,)^''™"' is (^^£±H£)1/p where A's are the largest eigenvalues 

m,,n— >oo 2g 

of corresponding T's, and the upper bound of lim /(m, rj,)^/™-" is {^2kY^'^'^ where (^'s are 
the largest eigenvalues of corresponding S's. Taking p = 2 , g = 4 and A; = 5, we get 

2.259132578... < lim /(m,^)^/™" < 2.259154406.... 

Note that the entropy constant of the generalized aztec diamond lattice is 

rj = lim f{m, 72)^/(2™"-+'"+"-) = lim f{m, where 2mn + m + n is the number 

m,n—*oo m,n— >cx) 

of vertices of vertices of generalized aztec diamond. Then we see that the entropy con- 
stant of the generalized aztec diamond lattice is ( lim /(m, 77,) i/'"") 1/2^ which is between 
1.503041110... and 1.503048371.... 



B32 = 



Thus 



^3 = 



4. Lattices with transfer multiplicity three 

The 8.8.4 lattice graphs Gm,n, as shown in FigJU are finite subgraphs of the 8.8.4 tiling 
of Euclidean plane. Some study of the properties of 8.8.4 lattice graphs can be found in 

Oils]. 



14 ZUHE ZHANG 




Figure 4. 



The cylindrical 8.8.4 lattice graphs are obtained by identifying the top row and the 
bottom row of the 8.8.4 lattice graphs Gm.,n- The toroidal 8.8.4 lattice graphs can be 
obtained by identifying the left cycle and the right cycle of the cylindrical 8.8.4 lattice 
graphs. 

Consider the transfer matrix of the 8.8.4 lattice graph Gm,n- Define the assembling of 
the {i + l)th column ViJ^i to the ith column Vi as step i as we did before. We can see that 
the generation of each independent set of Gm,n involves 3n — 3 steps after the first column 
Vi determined. Step one is to assemble V2 to the right side of Vi. The transfer matrix 
representing step one, denoted T^^ , can be constructed as follows. Let Rm denote the set 
of all possible vectors which can appear as Vi. Clearly consists of (2m — 2)-vectors of 
O's and I's in which no consecutive I's occupy the positions of the 2k — 1-th and 2k-th 
entries, for 1 < k < m — 1. Since there are three possibilities at each pair of consecutive 
2k — 1-th and 2A;-th 1 positions, the set Rm has 3*""^ vectors. The set of all possible 
vectors V2, denoted as Km, consists of all m-vectors of O's and I's. So Km has 2*" vectors. 
Thus the transfer matrix = {Tvi,v2) is a 3™~^ x 2™ matrix whose rows are indexed by 
vectors of Rm and columns are indexed by vectors of Km, where Ty-^^y^ = 1 if vi and V2 
represent possible consecutive pair of columns in an independent set of Gm,n and = 
otherwise. Similarly, the transfer matrix for step three is a 2™ x 3™"^ matrix T,„3 which 
is the transpose of T^^, and the transfer matrix for step two is a 2"^ x 2*" matrix 
whose rows and columns are indexed by vectors in Km, T^j's entry in position (a,/?) is 
1 if the vectors represented by a, (3 are orthogonal, and is otherwise. Note that T^^ is 
the transfer matrix for every step i when z = 3/c + l(0<A;<n — 2), is the transfer 
matrix for every step i when i = 3k + 2(0 < k < n — 2) and Tm-^ is the transfer matrix 
for every step i when i = 3A; (0 < A; < n — 1). Thus, the transfer multiplicity of Gm,n is 
3. Thus, if we take the transfer matrix of Gm,n to be T3=T3^T32T33, then it is a 2™ x 2™ 
symmetric matrix with transfer multiplicity one. 
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When m — 2, if we index the rows and columns in increasing order in binary numbers, 
then the transfer matrix of Gm^n is a 3 x 3 matrix given by the product of T2^, and 



L23- 



T2, 



/l 
1 



1 1 l\ 

1 

1 y 



-22 



/l 1 1 l\ 
10 10 
110 
10 



'-23 



fill 

1 1 
1 1 
1 



Thus To 



/ 9 6 6 \ 
6 3 4 
6 4 3 



Consider the transfer matrix of i^n,m which is obtained from Gm,n by identifying its left 
column and right column. The transfer matrix , which represents every (3A; + l)-th 
(0 < A; < n — 2) step, can be defined as a 3""^ x 2"~^ matrix of O's and I's as follows. The 
rows of are indexed by vectors of i?„ and columns are indexed by vectors of Kn-i, 
and the entry of in position (a, 13) is 1 if a, f3 represent possible consecutive pair 
of rows in an independent set of Hn^m, and is otherwise. Let denote the transfer 
matrix that represents every 3A;-th (0 < /c < n — 1) step. It is no difficult to see that 
i?„3 is the transpose of Bm- The transfer matrix i?„2 that represents every {3k + 2)-th 
{1 < k < n — 2) step is a 2"~^ x 2"~^ matrix whose rows and columns are indexed by 
vectors of Kn-i- The entry of Bn^ in position (a,/?) is 1 if a,/? are orthogonal, and is 
otherwise. Thus if we take the transfer matrix of -f^n.m to be S„ = Bn^Bn^Bn^, then it is 
a 2"~^ X 2"~^ symmetric matrix with transfer multiplicity one. 

When n — 3, if we index the rows and columns in increasing order in binary numbers, 
then the transfer matrix of i^n,m is a 8 x 8 matrix given by product of , B^^ and S33 , 

where 



Thus 
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Note that Gm.n has the same symmetric transfer matrix in both horizontal and vertical 
directions. So the same approach in [4] can be taken here. Then we can easily get the 
following results. The lower bound of lim /(m,?!)^/™" is (^^£±2£)1/p where A's are the 

m,n— >oo 2<j 

largest eigenvalues of corresponding T's. And the upper bound of lim /(m,n)^/'"" is 

m,n^oo 

(?2fc)^^^^ where ^'s are the largest eigenvalues of corresponding S's. Letting p = 1 , g = 4 
and A; = 4, we have 
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4.631583395... < lim /(m, n)^/™" < 5.765456528.... 

Let /(m, n) denote the number of independent sets of the 8.8.4 lattice graph Gm,n- Since 
the number of vertices of Gm,n is 4mn — 2m — 2n, by the similar reason as for the aztec 
diamonds, we can see that the entropy constant of the 8.8.4 lattice is 

( lim /(m,n)i/'"")i/^ 

m,n— >oo 

which is between 1.467007628... and 1.549560101.... 

Remark 1. In this paper we show that for many statistical models, the entropy 
constant of a plane lattice is the same as the entropy constants of the corresponding 
cylindrical and toroidal lattices. But this phenomenon may disappear for some other 
models. As shown in |18], for the dimer problem the entropy constants of cylindrical and 
toroidal lattices are different. 

Remark 2. For many statistical models, using a single transfer matrix is not enough. 
So, the approach and the concept of transfer multiplicity introduced in this paper can be 
used to deal with more complex statistic models. 
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